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Abstract 

We consider effective operators describing Dark Matter (DM) interactions with Standard 
Model fermions. In the non-relativistic limit of the DM field, the operators can be 
organized according to their mass dimension and their velocity behaviour, i.e. whether 
they describe s- or p-wave annihilations. The analysis is carried out for self-conjugate 
DM (real scalar or Majorana fermion). In this case, the helicity suppression at work in 
the annihilation into fermions is lifted by electroweak bremsstrahlung. We construct and 
study all dimension-8 operators encoding such an effect. These results are of interest in 
indirect DM searches. 



1 Introduction 



Significant experimental activity is currently devoted to the search for Dark Matter (DM) by looking 
at the excesses in cosmic ray production through DM annihilations (or decays) in the galactic halo. 
Detailed predictions for the fluxes largely depend on the particle physics model for the DM, e.g. 
the mass, the annihilation channels etc. It is therefore desirable to describe DM annihilations and 
their products within a general and model-independent framework and Effective Field Theory (EFT) 
provides such a tool. 

Of course, the EFT is only applicable whenever there is a separation of scales between the 
process to describe (the annihilation of non-relativistic DM at a scale ~ Mdm) and the underlying 
microscopic physics of the interactions (at a scale A). This may not always be the for 
example, in the case of supersymmetry with a compressed spectrum. 

If we want to describe the annihilation of two non-relativistic DM particles, whose relative 
velocity is v ~ 10~ 3 (in units of c) in our Galaxy today, it is convenient to expand the cross section 
in powers of v 

va = a + bv 2 + 0(t> 4 ) , (1.1) 

where the first term corresponds to annihilation in the state of orbital angular momentum L = 
(s-wave) while the second term describes L = 1 (p-wave). For the annihilation DM DM — > ff of a 
self-conjugate DM particle (real scalar or Major ana fermion) into SM fermions of mass m/, helicity 
arguments lead to a oc (mj/Mdm) 2 , and hence a very suppressed s-wave term for light final state 
fermions (e.g. leptons), while the p-wave is suppressed by v 2 . 

It is clear that the correct operator expansion must be done in terms of two parameters: mass 
dimension of the operator and relative velocity. The effective lagrangian would be generically given 
by an infinite series of non-renormalizable operators 

= Yl J^A (c (d ' s) (d ' s) + c (d ' p) (c! ' p) ) , (1.2) 

<i>4 

where cK rf ' sorp ) indicates that the operator of dimension d describes s-wave or p-wave annihilations, 
and we neglect annihilations in waves higher than p. There is no obvious ordering of the importance 
of the operators. An important role in this respect is played by ElectroWeak (EW) bremsstrahlung, 
which has received significant attention recently [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11] (for earlier studies 
on the impact of gauge boson radiation on DM annihilations or cosmic ray physics, see [12, 13, 14]). 
Taking into account processes with the inclusion of EW radiation eludes the helicity suppression 
and opens up an s-wave contribution to the cross section [2, 3, 6, 7, 10, 12]. 

In this paper we will classify the operators (up to dimension 8) according to the -u-behaviour of 
the amplitude connecting two self-conjugate particles - real scalars or Majorana fermions - in the 
initial state with the final state of two massless fermions and possibly a gauge boson. 

The effects of lifting the helicity suppression by means of EW radiation has so far been studied 
within the context of explicit models [2, 3, 6, 7, 10]. The results of the present paper provide a model- 
independent approach to this problem and can be used to place robust constraints on the new physics 
responsible for the DM sector. Although several analyses in the literature place phenomenological 
constraints on the coefficients of the dimension-six operators (see e.g. Refs. [15, 16, 17, 18]), the 
important role of higher-dimensional operators is typically underestimated. 
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The remainder of the paper is organized as follows. In Section 2, we explain our methodology 
and construct the effective operators contributing to s-wave DM annihilations. We compute the 
differential and total annihilation cross sections for each operator in Section 3, and compare them to 
the contribution from the typical lowest-dimensional operator. We conclude in Section 4, mentioning 
phenomenological applications and possible analyses that could be carried out using these results. 
Finally, we collect some useful relations and identities in the Appendices, for the convenience of the 
reader. 

2 Effective operators 

We will carry out our analysis under the following set of assumptions on the DM sector: 

1. the DM is either a Majorana fermion or a real scalar field; 

2. the DM is neutral under the SM gauge group; 

3. there exists a Zi symmetry under which the DM is odd and the SM is even; 

4. the DM couples only to the fermions in the SM spectrum, which are assumed to be massless. 

Of course, any of these assumptions may not hold in reality and if this is the case our analysis requires 
modifications. Assumption 1 specifies the conditions under which the helicity suppression is effective, 
while assumption 2 is there to simplify the discussion, even though it is not strictly necessary and 
relaxing this assumption can also lead to interesting effects (see Refs. [7, 19]). Assumption 3 is 
commonly used in DM phenomenology to ensure the stability of the DM particle. 

More clarifications about assumption 4 are in order. Considering mj ^ would introduce another 
mass scale into the problem and would render the operator classification much less transparent. 
Our analysis is still valid in the regime mf/K <C v, which may not hold for the third-generation 
quarks. The other piece of information in assumption 4 is that the DM particle only couples to 
the fermion sector of the SM. This needs not to be true, of course. Allowing for DM interactions 
with the other SM particles, namely gauge bosons or the Higgs, the possibility of having additional 
operators contributing to DM annihilations in s-wave opens up. For instance for Majorana DM, x, 
coupling to the Higgs doublet H or a generic field-strength F^ v one can have {xi^X){H^ H) and 
(x7 5 x) F^F^, which are CP odd, and (x7 5 x) F^F^ , which is CP even. For a list of effective 
operators connecting DM to vector bosons and/or Higgs bosons see e.g. Refs. [20, 21]. We will not 
consider these possibilities here. 

We want to classify the operators of dimension d = 6, 7, 8 according to the f-behaviour of 
the amplitude connecting two DM particles in the initial state with two massless SM fermions / 
and a gauge boson in the final state. We look for operators which are hermitian, gauge invariant 
(under SU(2)l ® U(1)y) and giving a non-zero contribution to the amplitude for the annihilation 
of DM into two fermions and one gauge boson. We will further classify the operators according to 
their CP transformation properties (see Appendix A for the relevant transformation properties). 
The operators containing ifif give no contribution to the process under consideration, due to the 
Equation Of Motion (EOM). In order to ensure manifest gauge invariance, we first introduce the 
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following notation for the covariant derivative 



^/ = (t^ + igT^ + ig'YfB^f. 
ffip = f(y„-igT a W^-ig'Y f B^ 



(2.1) 



(2.2) 



where T a = cr a /2, with a a being the usual Pauli matrices, and the charge Qf of the fermion / is 
related to its hypercharge Yf by Qf = T 3 + Yf. Furthermore, the field-strength Wn V of the SM 
gauge fields is related to the covariant derivatives as usual 



In the following, we will separately deal with the cases where the DM is a Majorana fermion or a 
real scalar. For simplicity, we will consider only left-handed SM fermions fi, but the analysis can 
be applied to right-handed fermions straightforwardly. 

2.1 Majorana fermion DM 

Let us first suppose the DM particle is a Majorana fermion x- It is possible to build several operators 
containing two DM fields, two SM fermion fields and zero or one gauge bosons. They can be built 
in full generality requiring gauge invariance and hermiticity, and further classified according to their 
CP properties. The Majorana-flip properties and the chiralities of the SM fermions make several 
structures identically zero. The only two Majorana fermion bilinears that are non-vanishing in the 
limit v — > are X7 5 X an d X7^7 X ( see Appendix B). As we are interested in s- waves only, we will 
limit our analysis to these two bilinears for the Majorana fermions. For the pseudo-scalar bilinear 
X7 5 X no contractions with SM fermions can be built, as they would vanish either by chirality or by 
the EOM of fermions, so we are left with the axial-vector bilinear. 

The lowest-dimensional terms that can be written out of two x' s an d two f^s are of dimension 
6. At this level, there is only one non-vanishing operator satisfying all criteria: 



The [i = component of this operator could a priori give a ^-independent contribution to the 
scattering amplitude, but it actually vanishes because of the identity u^(p)vf(— p) = 0. Therefore 
this dimension-6 operator only contributes to the p-wave, as noted in Ref. [3], due to the helicity 
suppression, which cannot be removed by simply radiating a gauge boson from the external final 
state leg. In order to look for s-wave terms, we need to consider higher-dimensional operators with 
one EW gauge boson whose radiation in the annihilation process lifts the helicity suppression. 

At dimenion 7, there are no operators contributing to the s-wave cross section. In fact, all 
possible structures vanish either due to Majorana-flip properties or the chirality of the SM fermions, 
or because of the EOM of the /^s. A priori, the fi = component of the operator (x7 5< 9^x) [/l7 m /l] 
would give a ^-independent contribution to the scattering amplitude, but it vanishes again because 
of the identity u^(p)vf(— p) = 0. 

At the level of dimension 8, there are several structures that can be built requiring gauge invari- 
ance and hermiticity. They contain two x' s ( m the bilinear X7 /i 7 5 x)? t wo /l's an d two covariant 



'^^u]f = {igT a )w«j, 
7^, XI = 7(igT a )w» . 



(2.3) 



(2.4) 




(2.5) 
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Name 


Operator 


CP 


0M1 


(ryVx) [(/iX) 


7^ 




+ 


0M2 


*e/ii/pa (x7 5 7 M x) 




+ 


CM3 


ie^pa (x7 5 7 M x) 


(/X) 7 p (^7l) - (i/X) 7 p (^7l)] 


+ 


CM4 


i (X7 5 7 M X) 




) - (7z^; 






0M5 


i (x7 5 7' 1 x) 


7 l7 a 







Table 1: List of dimension-8 operators contributing to the s-wave cross section for the annihilation 
of Majorana DM into two fermions and a gauge boson. 



derivatives. It is possible to reduce the number of independent operators, contributing to the cross 
section for the process under consideration, by using EOM and the identities in Appendix C. In 
addition, some structures can be related to each other by terms (like (x7 5 7 M x) 9 2 [/l7^/])j which 
do not contribute to the s-wave annihilation into two fermions and a gauge boson. Therefore they 
contribute in exactly the same way to the amplitude for the process we are interested in. There 
remain only five independent operators of dimension 8 contributing to the s-wave annihilation of 
DM into two SM fermions and a gauge boson, listed in Table 1 together with their CP conjugation 
properties. We remain agnostic about the presence or absence of CP violation in the Dark Matter 
sector, which can possibly induce CP violation in the SM at loop level and therefore be further 
constrained. 

All other operators have either a larger dimensionality or produce more powers of v 2 in the 
annihilation cross section. Notice that we chose to keep a Lorentz-covariant formalism, despite 
looking at the non-relativistic limit. This implies that the same operator can lead to both v- 
independent and ^-dependent terms in the amplitudes; for example, the operator Omi also gives a 
contribution to the p-wave cross section, but we will not consider it as it is very suppressed. 

To summarize, the dimension-8 operator contributing to the s-wave annihilation cross section 
of Majorana DM into SM fermions is given by the sum of the operators in Table 1, 0™ = 
Y^i=i c i^®M.i\ the leading interactions are therefore described in terms of only a few operators 

1 1 5 

^eff = j2 c(6 ' P) °M + X4 ^ c i °™ + higher-dim . (2.6) 

i=i 

In absence of CP-violation in the DM sector, only the first three operators in the sum need to be 
considered. 

2.2 Real scalar DM 

Next, let us consider the case where the DM particle is a real scalar cf). By angular momentum 
conservation, two real scalars cannot annihilate into two massless fermions in the configuration with 
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Name 


Operator 


CP 


o m 








+ 




if 






+ 


OR3 








+ 


OR4 




+ 


0R5 


</> 2 [ 








Om 


*0 ^pvpa 








Om 


^0 £pVp<T ( 







Table 2: List of dimension-8 operators contributing to the s-wave cross section for the annihilation 
of real scalar DM into two fermions and a gauge boson. 



L = 0; this process would only occur through a chirality flip induced by a mass term. For this 
reason, in the limit mj = we are considering here, no s-wave annihilation cfxj) — > /l/l is possible. 
One can recover this result in the language of effective operators. At dimension 5, only a single 
operator can be constructed which, however, vanishes due to chirality, (f> 2 {f lIl) = 0. At dimension 
6, we have (p 2 d^(f lI^Il) = 0, by the EOM (cf. Eq. (C.2)). At dimension 7, there are no possible 
Lorentz contractions to construct an operator. 

Nevertheless, at the level of dimension 8, several gauge invariant hermitian operators can be built 
out of two c/>'s, two f^s and covariant derivatives. As discussed already for the Majorana case, it 
is possible to reduce the number of independent operators, contributing to the cross section for the 
process under consideration, by using the EOM and the identities in Appendix C. We are left with 
four CP-even operators and three CP-odd operators, listed in Table 2. A v-dependent annihilation 
(f)(f) — > fhfh is mediated by the operator Oru, while the s-wave annihilation of two fts can proceed 
by switching on the emission of a gauge boson in the final state. 

Other operators involving 4>d^(j) or a gauge field-strength can be obtained from the listed ones 
by integration by parts, using the EOM of the DM particle or the identities (2.3)-(2.4). For instance, 
the operator fd v [/l7 m /l] F^ u , considered in Ref. [22], is expressed in this basis as l/g(0Ri — Or2)- 

To summarize, the dimension-8 operator contributing to the s-wave annihilation cross section 
of real scalar DM into SM fermions is given by the sum of the operators in Table 2: 0^'^ = 
Sj=i Ci" S ^0~Ri] the leading interactions are therefore described in terms of only a few operators 

^fi = X4 E c T' S) °Ki + higher-dim . (2.7) 

i=l 

In absence of CP- violation in the DM sector, only the first four operators in the sum need to be 
considered. 
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,V(k) 





Figure 1: Diagrams for the annihilation process in Eq. (3.1). 



3 Annihilation cross sections 

In this section we show analytical results for the annihilation cross sections due to the operators 
found above. For simplicity, we restrict ourselves to considering left-handed SM fermions only, but 
the results can be easily adapted to account for annihilations into right-handed fermions as well. We 
consider the process 

DM(fci) DM(A; 2 ) -> f iiL ( Pl ) f j:L (p 2 ) V(k), (3.1) 

where light fermions in the final state - described here by the generic SU(2)l doublet F = (fx, f 2 ) T 
- can be both leptons and quarks, and where V = W^, Z,j. Note that diagrams with gauge boson 
emission from the final state legs have to be included in order to compute a gauge invariant amplitude 
(see Fig. 1). 

It is convenient to introduce the kinematical variables y, z defined by 

pS = (l-tfh/I/2, (3.2) 
p° 2 = (l-*)VS/2, (3.3) 
k° = (y + z)^/2, (3.4) 

which are subject to the following phase space constraints 

2 2 2 

mi, mf, mfr 

— <y<i, — < z <l-y+ — , (3.5) 
s sy s 

where s = (ki + k 2 ) 2 = 4M^ M /(1 — v 2 /4). The scattering amplitudes will be proportional to a 
coefficient A containing the correct gauge couplings according to the different possible final states 
and their SU(2)l (g> U(l)y quantum numbers; more explicitly, 

A(hfa) = -^f(i + y f ), A(f 2 f 2l ) = 

A(hhZ) = -JL-[l-(l+ yf y w ], A(f 2 f 2 Z) = +JL.[l-(l- yf ) l j* r ], (3.6) 

A(hf 2 W+) = -j=, A(f 2 fxW-) = -^=, 

where Yf = y//2 (e.g. yi = —1, yq = 1/3), g is the gauge coupling and s\v, cw are the sin and cos 
of the weak angle, respectively. 

The double-differential annihilation cross section is given in terms of the squared amplitude |.M| 2 
(averaged over the initial spins and summer over the final ones) as 

d 2 a \M\> 



dydz 128 vr 3 



(3.7) 
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while the total cross section is obtained by integrating over the kinematical variables on the phase 
space domain defined by Eq. (3.5). Let us now show the results of the computation of these cross 
sections for each of the operators found in the previous section. 



3.1 Majorana fermion DM 



In the zero-velocity approximation, s 
annihilation process (3.1) mediated by the dimension-8 operators Omi 



4Mq M and the double-differential cross sections for the 



, 0M5 in Table 1 are 



d 2 a 



dydz 
d 2 a 



dydz 
d 2 a 



dydz 

d 2 a 

i 

dydz 

d 2 a 



dydz 



0M3 



0M4 



„(8, S ) |2 ^M 
'Ml I 27r 3 A 8 



1 - y - z + 



1117 



4< M 



2 i 2 

y +z 



in 



v 



?M 2 



„(M|2 
"M2 I 

JM|2 
-Ml I 

„(M|2 



"Ml 



d 2 a 



dydz 
d 2 a 



dydz 



Omi 



-M4 



7T 



3 A 8 



z)(y 2 + z 2 \ 



DM 



2(y + z) + 2] 



d 2 a 



„(M|2 

-M5 I 

„(M|2 U dydz 
-M4 I 



(3.8) 
(3.9) 
(3.10) 

l" > 
(3.11) 

(3.12) 



Notice that in the limit my/MoM —> the differential cross sections are the same for all operators, 

up to an overall numerical factor. 

In the indirect searches for DM, the observables measured experimentally are the fluxes of cosmic 

rays, which are directly related to the energy spectra of particles generated by DM annihilations 

at the production point. By integrating the double-differential cross sections listed above once, one 

obtains the energy spectra of the SM fermions and of the gauge bosons at production. We consider 

the distributions of the final fermion energy (EA and of the final gauge boson energy (Ey), defined 

as _ 

dN 1 At(DMDM — > ffV) 

dlnx cj(DMDM^//y) dlnx 
where x = Ejy /Mdm, and shown in Fig. 2. 

The distributions originating from the set of operators Omi, 0M2> Om3 differ just by an overall 
factor, which cancels out by normalizing the spectra, and therefore produce the same curves. The 
same argument applies to the other set 0m4,Om5- However, the operator Omi would instead give a 
different result with respect to Om4 ; but the difference is not visible in the fermion energy distribution 
(left panel of Fig. 2) because of the smallness of my /Mdm- 

The dimension-6 operator Cm mediates the two-body annihilation \X ~^ IlIl m p-wave. The 
processes where a gauge boson is radiated from the final state fermion are still in p-wave, as dis- 
cussed in Ref. [3]. Therefore, the energy spectra of fermions and gauge bosons originating from 
the dimension-6 operator are going to be much less important than those from the dimension-8 
operators. The total cross section for the two-body process XX ~^ JlIl mediated by Cm is simply 



M 2 

I J6,p)|2 DM 2 

1 1 12vrA 4 



(3.14) 
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■S3 



Ml, M2, M3, M4, M5 



X = £ f /M DM 



M DM = lTeV E 




X = E V /M DM 



M DM = 1 TeV 



Figure 2: Majorana fermion DM. The energy distributions dN j dln(E / -Mdm) of the final fermion 
(left panel) and of the final gauge boson (right panel), for the different dimenion-8 operators. We 
set = c^ s) = 1 and M DM = 1 TeV. 



while the total cross sections of the three-body processes (3.1) mediated by the dimension-8 operators 
are obtained by integrating over the full phase space 



va 



va 



va 



va 



va 



Ml 



M2 



M3 



Mi 



M5 



-Ml 



s),2 A Mp M 
240vr 3 A 8 

„(M|2 
-M2 I I 
„(8,.) |2 W lMl > 
-Ml I 
M)|2 



15 



I17T 



Km 



+ 5 



ml 



Km 



. 4 + 6lQ 2 Mm) +0 

my ) 



777 



V 



Km 



-M3 



„(8,a), 2 



va 



Ml i 



I°M1 

' M4 1 120vr 3 A 8 

„(«.«) |2 



8 + 15 



mf 



^DM 



+ 



V 



Km 



40 



my / 



<M 



-M5 



J8,s)| 



Ml 



,(3.15) 
(3.16) 

(3.17) 
,(3.18) 
(3.19) 



-M4 



For simplicity, we have only reported here the leading terms in the expansion in powers of my/MoM, 
but the complete analytical expressions are used in the plots. The sub-leading terms can be of the 
same order as the contributions from higher-dimensional operators we are neglecting. 

The relative importance of the s-wave three-body process due to dimension-8 operators with 
respect to the p-wave two-body annihilation due to the dimenion-6 operator is captured by the 
ratios of the total cross sections, plotted in Fig. 3. Three-body cross sections can be sizeably larger 
than the two-body ones. It is evident that the dimension-8 operators dominate the total cross section, 
provided that the effective operator scale A is not too large with respect to the DM mass Mdm- It 
is clear that limiting an EFT analysis for DM annihilations to the dimension-six operators misses 
the right result, as the cross section receives important contributions from operators of dimension 
higher than six. 

One could have expected this result by an order-of-magnitude estimate of the two-body and 
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500 1000 1500 2000 2500 3000 2 4 6 8 10 



M DM [GeV] A/M DM 

Figure 3: Left panel: The ratio of the total cross sections va(xx ~~ ^l^L ^) Mi/ V<J (XX ~~ ^l^l)|m 
as a function of Mdm- We set c^ 6,p ^ = c[ S ' s ^ = 1, and A = 5 TeV. The curves corresponding to 
Cm4jC'm5 (i n green) departs from those o/0m2;0M3 at low M-qm and are slightly higher. Right 
panel: The ratio of the total cross sections vcr(xx ^~l^l^)\mi I V<J (XX ~^ ^l^l)Im as a f unc ti° n 
o/A/Mdm- We set c^ = cf' s) = 1 and M DM = 1 TeV. 



three-body total cross sections originating from the operators above 

va(2^2)\ 0e ~ |c^|V^M (3.20) 

va(2^3)\ 0a ~ I^.JI——-^-, (3.21) 

where ajy = g 2 / (47r), from which we learn that the three-body cross section due to the dimension-8 
operator can be bigger than the two-body one due to dimension-6, provided that 

A 2 „ 1 fa^ 



| c (M|m Sm ^W^" S ° (3 - 22) 

for t> = 10~ 3 . This estimate is confirmed by the numerical results in Fig. 3. 

Let us conclude this subsection with a comment on the toy- model studied in Ref. [3], consisting 
of a Majorana spinor x with mass Mdm and a scalar S with mass Ms > Mdm, in addition to the 
SM particle content. The added particles are a SM singlet and a 5C/(2)i-doublet respectively. The 
Lagrangian is of the form 

C = C SM + \x(i$ ~ M DM ) X + (D^S)\D^S) - M%S*S + (y x x(Fia 2 S) + h.c.) . (3.23) 

In the limit Ms S> M^Mj the interactions can be accurately described by the effective operator Omi 
in Table 1. Indeed, in the limit v — > 0, the amplitudes for the process XX ~^ IlIlZ as due to the 
toy model Eq. (3.23) and the effective operator Omi match, provided that 4cjJq /A 4 = —y^/Mg. 

3.2 Real scalar DM 

As discussed in the previous section, there is no two-body annihilation <p(f) — > fifi m s-wave in 
the limit mj = 0. So the EW bremsstrahlung opens up a three-body annihilation channel which 
is otherwise absent. The first non- vanishing contribution to the s-wave annihilation cross section 
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comes from the dimension-8 operators in Table 2 and there is no contribution from lower-dimensional 
operators to compare with. 

For the case of real scalar DM, we computed the double-differential cross sections for the anni- 
hilation process (3.1) mediated by the dimension-8 operators Ori, . . . ,Or7 listed in Table 2 



d 2 a 



dydz 



d 2 a 



JM|2 4,4 2 M DM 



dydz 



d 2 a 



Rl 



R2 



-Rl 



dydz 



d 2 a 



dydz 
d 2 a 



dydz 
d 2 a 



dydz 
d 2 a 



dydz 



R3 



R4 



R5 



R6 



R7 



„(M|2 
-R2 I 



„(M|2 
-R3 I 



7T 3 A 8 



7T 3 A 8 



A 2 M 8 
A DM 

TT 3 A 8 



;i - y - z)(y 2 + z 2 ) + 



(1 - y - z){y 2 + z 2 ) 



ml 



[(y + z) 2 - 2{y + z) + 2] 



ml 



4< M 



[3(i/ + z) 2 -2{y + z + 3yz) - 2] 



ml 



+ 



ml 



16M DM 



(3.24) 



(3.25) 



l-y-z+ 



m 



v 



*M 2 M 



2 i 2 

y +z 



ml 



2M 2 
DM 



[2(y + z)-l} + 



ml 



4M DM 



1 c 



Rl 



16|J 8 ' S )|2 dwdz 
l c Ri I y 

(M 2 A 2 M 8 M 

C R5 I 

\M\2 
l L R6 I 



Rl 



7T 3 A 8 

d 2 a 



l-y-z+ 



in 



v 



4M 2 M 



2 i 2 

y +z 



ml 



9 /If 2 



-R5 



; > s )|2 



-R7 



s«)|2 



dydz 
d 2 a 



( 8 > s )|2 dydz 



-R5 



R5 



K5 



(3.26) 
(3.27) 

(3.28) 

(3.29) 

(3.30) 



They are clearly symmetric under the exchange y o z. The energy distributions of the fermion and 
gauge boson, defined as in Eq. (3.13), are shown in Fig. 4. Some of the distributions (e.g. those 
due to the operators Ori, Or4) differ just by an overall factor, which cancels out by normalizing the 
spectra, and therefore produce the same curves. In the limit my <C Mdm the distributions from 
all different operators become proportional to each other. Therefore, because of the smallness of 
my /Mdm, the differences in the fermion energy distribution are not visible (left panel of Fig. 4). 

For convenience, we also report the results for the total cross sections, obtained after integrating 
over the kinematical variables 
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Figure 4: Real scalar DM. The energy distributions dN / dln(E / M~e>m) of the final fermion (left 
panel) and of the final gauge boson (right panel), for the different dimenion-8 operators. We set 
c^; s) = 1 and M DM = 1 TeV. 
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,(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 

(3.37) 



As for the Majorana case, we have only shown the leading terms in the expansion in my/McM- The 
sub-leading terms can be of the same order as the contributions from higher- dimensional operators 
we are neglecting. 



4 Conclusions 

We have considered the annihilation of a self-conjugate DM particle - either a Majorana fermion or 
a real scalar - into light SM fermions. In these cases, the simple annihilation DM DM —> Jl/l 1S 
helicity suppressed. The radiation of electroweak gauge bosons lifts the suppression and opens up 
an s-wave channel, independently of the relative velocity of the annihilating DM particles. From the 
effective operator point of view, we found that this effect is encoded by dimension-8 operators. 



11 



The dimension-8 operators, despite suffering from a high mass dimensionality, enjoy no sup- 
pression by the relative velocity. We have found all the dimension-8 operators mediating s-wave 
annihilations of DM into SM fermions and a gauge boson, and for each of them we have com- 
puted the differential and total cross sections. These operators encode the so-called "virtual internal 
bremsstrahhmg" effect. 

In the Majorana DM case, the two-body annihilation DM DM — > /l/l proceeds through a 
dimension-6 operator and is velocity suppressed (p-wave). We have shown explicitly and quantita- 
tively that the dimension-8 operators actually provide a bigger contribution to the annihilation cross 
section than the dimension-6 operators, for DM masses well above the weak scale. 

For real scalar DM, the s-wave annihilation into massless fermions is forbidden, and becomes 
allowed in the presence of an extra gauge boson in the final state. We found the dimension-8 
operators mediating annihilations into SM fermions and a gauge boson, which switches on an s- 
wave annihilation process. 

In conclusion, our analysis has strengthened the statement that, when dealing with the annihi- 
lation of heavy self-conjugate DM today, it would be erroneous to neglect the effects of dimension-8 
operators. The results presented in this paper are of interest also when applied to phenomenological 
analyses, as they provide a model-independent language into which to translate the experimental 
data, present exclusions and possible future signals. 
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A CP transformation properties 



Under charge conjugation the spinor fields if) and ip transform as 



(A.l) 
(A.2) 



where the charge conjugation matrix C obeys C T = = C _1 = — C. Under parity transformation 



if) -> P~VP = i/3i/j , 
ip ->• P-^P = -itp/3 . 

For instance, in the Dirac representation of the 7-matrices one has 



(A.3) 
(A.4) 



C = i7 2 7° 



A 



-10 



-10 



/3 = 7° 



1 
-1 



(A.5) 
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The following identities turn out to be useful 

7°7 2 7^L7 2 7° 

where 

+1, a* = 0, 



(-iy 



(A.6) 
(A.7) 



(A.8) 



-1, M = l,2,3. 



Using these relations it is easy to prove the following transformation properties under CP of operators 
built out of left-handed fermion fields fi ■ 





cp{ 


'h-f^fL + fL^-ffL 




(-i)(-in-ir, 


(A.9) 




CP\i 






(-i)^(-i)^ 


(A.10) 




cp{\f L H> 




} - 




(A.11) 


CP {i [fj) (t p f L ) - (/ L ^ M ) %f L 


} - 


(-i)m. 


(A.12) 


CP{ 




} - 




(A.13) 


CP\i 


f L J^ p (t<>f L ) - (j L fi p ) fr>ff L 






(A. 14) 






CP{f Llv W^f L ] = 


(-I)a*. 


(A.15) 






CP{j Llv W^f L ) = 




(A.16) 



If X is a Majorana fermion field, the bilinears XJ 5 X: X7 M 7 5 X transform under CP as 

-1, 

while the bilinears built out of a real scalar <fi transform as 

i-in-iy. 



cp{xi b x] 
cp{x^i 5 x} 



CP{<t>d^} 
CP{d^d v 4>] 



(A.17) 
(A.18) 



(A.19) 
(A.20) 



B Majorana fermion bilinears 

The standard decomposition of a Majorana fermion field in momentum space is 

X (x) = [ - f * Y, U(k)u s (k)e-^ + at(*M*)e* x ) , (B.l) 
J (2vrj \/2h k g \ / 

from which one can compute the matrix element between the initial state with two DM particles 
and the vacuum, for the bilinears at small velocities. In the Dirac representation of 7-matrices we 
obtain 

<0| xO X \X(P°,P)X(P°, ~P)) ~ v(-p)Ou(p) - v(p)Ou(-p) , (B.2) 



13 



where the Majorana nature of the initial particles implies to consider the process with spins and 
momenta interchanged (and a relative minus sign due to Fermi statistics). Let us first consider the 
case where the operator O is one of the basis matrices Ta = {1, 7 5 , 7^, 7 5 7^, a^ u }. Noting their 
behaviour under charge conjugation 



(r A = M 7 5 , 7 V) 
(r A = 7 ",o 



(B.3) 



we find that 



v(p)T A u(-p) = < 



-v(-p)T A u(p) (T A = 1, i 7 5 , 7V) 
+v(-p)T A u(p) (F A = ^,a^) 



(B.4) 



hence the vector and tensor operators do not contribute to the annihilation of Majorana particles 
- as it is well known - and it is sufficient to compute the bilinears for scalar, pseudo-scalar and 
pseudo- vector interactions 



v(—p)u(p) = —2rf{a-p)^ ~ v 

v(-p)"/°j 5 u(p) = (E + M) »jt£ 







~ V 







v(-p)f-/ 5 u(p) 



The bilinears u(p)T A v(— p) can be obtained from the previous ones by 



u(p)T A v(-p) = v(-p)(l°T A J>(p) 



such that, for example 



u{p)l°l 5 v(-p) = (E + M)^ V , 
u(p)7 i 7 5v (-P) = 2i(pxaY^r]. 

If the operator O contains one derivative, we have 

XJ 5 d x ~ v°, 
Xl^diX ~ v, 
X7 5 7 ^oX 

% 



xr'fdjx = 0. 



(B.5) 



(B.6) 
(B.7) 



(B.8) 



The generalization to an arbitrary number of derivatives is straightforward. 
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C Useful Identities 

By a direct computation we find the following identity involving one covariant derivative 

du QlI^l) = 7i7 M + Q L %)-ffL • (ci) 
A particular case is the contraction p = u, which gives, using the equations of motion 

d, (7l7 m /l) = . (C.2) 

For two covariant derivatives we find 

d p d v (f L -ff L ) = J L r^ P 0uf L ) + (jfiph^ufL) + (7lK)7 m A/l) + (fL^ufip-fh ■ 

(C.3) 

A useful relation can be extracted by contracting p = p 

(7l7 m /l) = fL$0ufL) + (7lK)S/l • (C.4) 
However, exchanging the partial derivatives on the left hand side, and using Eq. (C.2) we find 

J L ~$0 v f L ) + (f L K)*0f L = . (C.5) 

Eq. (C.3) shows that the action of a partial derivative is equivalent to the action of the corresponding 
covariant derivative on each single term. This relation is crucial to relate operators containing 
quadratic (d 4> 2 ) and derivative (D (pd^fi) terms of the scalar field. Notice that this kind of relation 
remains true also for differences of SM operators instead of sums. By direct computation we find 



dp 



J Ll ^ p v f L ) + (f L t p )-f(B u f L ) - (7lK)7 m P h) - (fifrvjbp-ffL ■ (C6) 
Similar relations hold for three covariant derivatives 



d a d p a v (7l7 m /l) = J L ^t a t p u f L ) + (f L ^ayft p v f L ) + Q L t p yfi a v f L) 

+ (f L t p )^ a ^0 u f L ) + (7lK)7 m ^(^ p /l) + (7l^,)K7 m (^p/lJ 
+ (7 L K)^p7 M (^/L) + (7LK)^ P K7 At /L • (C7) 



This cumbersome identity simplifies contracting p = p,; using the same argument previously exploited 
we find 



+ (f L fivY0'f(l5*fL) + (7l^,)5K7 m /l , (C.8) 



which is useful to find dependencies among seemingly independent operators. 
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